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Abstract. We prove all the maximizers of the sharp Hardy-Littlewood-Sobolev 
inequality are smooth. More generally, we show all the nonnegative critical 
functions are smooth, radial with respect to some points and strictly decreas- 
ing in the radial direction. In particular, we resolve all the cases left open by 
previous works of Chen, Li and Ou on the corresponding integral systems. 



1. Introduction 



The classical Hardy-Littlewood-Sobolev inequality states that for < a < n, 
hat -i- + - 

Po 5 

f(x)g(y) 



1 < po, qa < ^ such that ^ + ^ = 1 + f (see [S] theorem 1 on pll9]) 



/ 



\x - y\ 



-dxdy 



<c(n,p ,a) |/| LP0 (R") \9\ L "omn) ■ 



In [EI], it was shown that the sharp constant 
c{n,p ,a) = sup 



a dxdy : |/| iP0 ( Rn ) 



1. Ifll 



f - y\ 



L90(K") 
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is achieved by some functions / and g. Moreover, after multiplying some constants, 
any maximizer /, g must be radial symmetric with respect to the same point, strictly 
decreasing in the radial direction and satisfy the integral system 



/(*) 



Po-l 



g(y) 



f - y\ 



9o-l 



f{y) 



It was also shown that when po = 9o> we have 
f[x) =g(x) = c(n,p ) 



X 



f - y\ 



Vpo 



dy. 

n—a & 



\x — x \ + A 2 



for some A > and xq S R™. 



If we let p = q 
equation becomes 

(1.1) u(x) = 



u — f p " 1 , v = g q ° 1 , then the Euler-Lagrange 



v{y) q , , , 
— -dy, v (s) = 



f - y\ 



f - y\ 



and v e Li +1 (R n ) and < a < 



for nonnegative functions it € L p+1 ( 

— — < p, q < oo, -4t H — rr + - = 1- When p = q = 2±« as observed in O, it 

n—a 5 p+1 q+1 n * ^ n—a 1 & 

follows from the fact ,*_„ = c (n, a) — ^— s- * — ^— s- that u = v, then the system 
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reduces to 



(1.2) 



u 




dy. 



In |CLQ2| . using an integral form of the method of moving planes ([GNN]), it was 
shown that any nonzero nonnegative regular solution u of (|1.2|) must be of the form 



for some A > and xq £ K". This solves an open problem proposed in [Li] (see a 
somewhat different argument in [L] and the clarifications in [CLQ21 Remark 1.3 on 
p332]). In |CL[ ICLOT] , such kinds of analysis were extended to the system (II. ip 
under the additional constraints p > 1 and q > 1. However the analysis does not 
give the regularity of maximizer for all the Hardy-Littlewood-Sobolev inequalities. 
On the other hand, it does not seem that we will have nonsmooth maximizers for 
the Hardy-Littlewood-Sobolev inequality in any case. The main aim of this article is 
to prove the regularity and radial symmetry of nonnegative solutions of the system 
(|l.ip in its full range. Another motivation comes from the study of regularity issues 
for a similar integral system in |HWY], 

Theorem 1.1. Assume < a < n, — — < p,q < oo, 4r H rr + — = 1. 

u £ L p+l (M™) is nonnegative and does not vanish identically. If 



Then u £ C°° (M n ), v £ C°° (R n ). Moreover, there exists a point x £ K" such that 
both u and v are radial symmetric with respect to Xq and strictly decreasing along 
radial direction. 

Indeed, the regularity is still true under the relatively weaker assumption u £ 
Lf+^M") (see Proposition . The method in [CTI ICLOlj . which is basically 
linear in nature, does not seem to work for the case when one of the two indices 
p and q is strictly less than 1. We will develop some nonlinear approaches which 
work for all p and q at once. In HWY] , we will apply this technique to derive the 
regularity for another integral system. In Section [2] below, we will prove a local 
regularity result which has the regularity part in Theorem 11.11 as a corollary. In 
Section [3) we will prove all the solutions are radial. 
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In this section, we will show any solution u, v to the system (|1.1[) must be smooth 
if we assume u £ L^ 1 (R ra ). Such a local integrability condition is necessary for 
the smoothness because as observed in [Li], system (jl.lj) has singular solutions as 





2. Regularity issue 



u(x) — c (n, a,p) \x\ 



p+i 



v (x) = c (n, a,p) \x\ q+1 . 
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This follows from a simple change of variable in the integrals. To achieve the 
regularity, we start with a local result which has some similarity to [CJLI4 theorem 
2] and p2 theorem 1.3]. 

Proposition 2.1. Given < a, (3 < n, 1 < a, b < 00, 1 < r < 00 such that 

1 1 _ a (3 
ra b rn n 

Assume 

n 

< p < q < 00, 

n — p 

a r 1 r 1 
-<- + -<- + -< 1, 

n q a p a 

u, f G L P (B B ), U G L a {BR), V G L b (Br) are all nonnegative functions with 
f\ BR/2 EL«(B R/2 ), 



TO 

and 

l/r 



U (a;) < 



V(y) 



b r \x - y\ 



I in— a 

b r \y - z\ 



dy + f (x) 



for x G Br, then u G L q (_Bfj/ 4 ) ; moreover 

\u\ L v( Br/4 ) ^ c{n,p,q,r,a,[3,a,b) (R*~9 \u\ Lp(Br) + I/I l ,( Br/2 )) • 

Proof. By scaling, we may assume R = 1. First assume we have u, f E L q (B\). 
Denote 



Jb, \x-y\ 



Ib 1 \x-y\ r 
Let pi and (?i be the numbers defined by 

1 r 1 a 1 r 1 a 
Pi p a n qi q a n' 
then it follows from Hardy-Littlewood-Sobolcv inequality that 

\v\ LP1(Bl) < c(n,p,r,a,a) \U\ La(Bi) \u\ LP(Bi] , 

\v\ L n( Bl) < c(n,q,r,a,a) \U\ La(Bi) \u\ r Lq(Bi) . 

Given < s < t < 1/2. For 1 6 B s , we have 

»<■'■' / ' - , •"''/ / V , iy)v{ ^ r dy + f(x) 

b 1 \b 1±± \x-y\ 

2 

< / '."""HI. *v + / f fe) - (»)"' * + / w 



V(y) 


v (y) 1/r 


\x- 


y\ n - p 


V(y) 


v (y) 1/r 


\x- 


y\ n - p 


V(y) 


v (y) 1/r 


\x- 


y\ n - p 



Hence we have 

l/ r C ( n )Pl <7i r : a i Pi fl j Mz,p(B ) 

< c(n,q,r,(3,b) \V\ Lb(Bi) M + x+ _ n _ -+l/lz,«(s 1/2 ) 
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On the other hand, for x G B s+t , we have 

2 

U{y)u(y) r 



v {x) 



r u( y )u( y y dy r 

JB t \x-y\ n ~ a JB 1 \B t 
U(y)u(y) r c(n,a) 

I in — a y ' 

B t \x-y\ 



U (y)u(y) r 

I in — a y 

f - y\ 



< 



„ „ , U(y)u(y) r dy 
(s - t) JB 1 \B t 

U (y) u{y) r c(n,p, r, a, a) |f/| L a (Sl) |u|£, (Bl) 



This implies 



< c (n, g, r, a, a) \U\ La(Bi) \u\ r Lq(Bt) + 



Combine the two inequalities together, we see 



/ ,\n — a 
(s - t) 

c(n,p,q,r,a,a) \U\ L « {Bl) \u\ r LP{Bl) 
(s-t) 



\Li{B 



s) < c(n,q,r,a,/3,a,b) \U\ L ^ {Bi) \V\ Lb{Bi) \u\ Lq[Bt) 
c(n,p,q,r,a,P,a,b) 

( s _ f \max{(n-a)/r,n-/3} l U li p (Bi) + !■> li« (B1/2) 

< 1 I I c(n,p,q,r,a,f3,a,b) 

- 2 |U| i"(S t ) , _ ^max{(ra-a)/r,n-,S} l W liP(Bi) + 1/ lid (B 1/2 ) 

if e is small enough. It follows from usual iteration procedure ([HL, lemma 4.3 on 
p. 75]) that 



Mi«(B 1/4 ) <c{n,p,q,r,a,f3,a,b) (\u\ LP[Bl) + \f\ L ^ Bl/2 ) / 
To prove the full proposition, we note that for some function < rj (x) < 1, 



u (x) = r\ (x) 
We may define a map T by 
T ((p) (x) = 77 (x) 
Note that we have 



V(y) 



x-y\ 



u(z)u( z y 



r \y~ z \ 



dz 



l/r 



dy + T)(x)f (x) 



V(y) 



Br \x - y\ 



n-0 



1 in — a 

b r \y-z\ 



dz 



l/r 



dy. 



I t (v)Ilp(Bi) ^ c(n,p,r,a,p,a,b) \U\^ Bi) \V\ L „ {Bi) W\ LP(Bl) < ~ \v\ LP{Bl) 



and 



\T{f)\ Lq(Bl) <c(n,q,r,a,0,a,b) \U\]^a {Bi) \V\ L „ (Bi) \<f\ Lq{Bl) < | \<p\ L «( Bl ) 

if e is small enough. Moreover, for any ip,ip £ LP (Bi), it follows from Minkowski's 
inequality that 



2 \^\LP(B! 



1 



hence 



|T(vj)(a:)-T(^)(x)|<r(|vJ-V|)(a;) fcrie^, 



i p (-Bi) 
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Similarly, we have for any ip,ip £ L q (i?i), 

\T {(p) - T W\ Lq{Bl) <\W- *I>\l<(Bi) ■ 

For k £ N, let fk (x) = min{/ (x) , fc}, then it follows from contraction mapping 
theorem that we may find a unique Uk £ (-Bi) such that 



u fc (x) = T (u k ) (x) + 77 (x) / fe (x) 



7?(x) 



V(y) r ' rr^i...^ir 1 Vr 



1 \ n ~ 
b r |x - y\ 



U{z)\u k {z)\ r 

1 in — a 

|y - 2I 



o!y + 7? (x) / fe (x) . 



Applying the apriori estimate to it/., we see 

|wfcl £ «(B 1/4 ) <c(n,p,g,r,a,^,a,6) (\u k \ LP{Bl) + \f\ L i(B 1/2 )) ■ 
Now observe that 

u (x) = T (u) (x) + 77 (x) / (x) . 

We see 

\u k - u\ LP[Bi) < \T{u k ) - T(u)\ LV[Bi) + \ f k - f\ LP(Bl) 

< g K ~ W liP(fl!) + l/fe - /Ilp(S!) ' 

Hence \uk — u\ LP ^ Bi ^ < 2 — /Ilp(b 1 ) ~~ * as fc — » 00. Taking a limit process in 
the apriori estimate for Uk, we get the proposition. □ 

Now we are ready to derive the full regularity for the system (jTTTJ) . Such regu- 
larity under the additional assumption p > 1 and q > 1 was proved in |CL[ |L] . 

Proposition 2.2. Assume < a < n, —2— < p,q < 00, — tt H rr + — = 1, 

it £ ifo^ 1 (R n ) is nonnegative and does not vanish identically. If 

.<•>-/ r 5 ^* •«-/ r^*' 

je" |x - y| jk« |x - y| 

TTien u,u £ C°° (R n ). Moreover i/ we hora u £ (R n ), rj/ien u (x) -> and 

v (x) — > as |x| — > 00. 

Proof. Since u £ (R"), we see it(x) < 00 a.e. x £ R n . It follows that 

w (x) < 00 a.e. x £ M™. For any i? > 0, we may find xo £ Br such that v (xq) < 00. 
This gives us L„, „ 1 — ofy < 00. It follows that L„> „ rr^-^dy < 00. Now 



1 \ f u fa) ^ , / M (y) . 

JSi? F - y| JR"\s fl |x - y| 

it follows from the Hardy-Littlewood-Sobolev inequality that the first term lies in 
L q+1 (R n ). On the other hand, for x £ Bqr with < 6 < 1, we have 



/ 

■J 7 j 



u ( yf , < 1 /" u (yf 

\n — a o — /-x n \n—a I i in — a "* 



\b b \x - y\ ' (1 - 6) Jr»\b r \y 



It follows that v £ L q ^} {Br). Since R is arbitrary, we have v £ L q ^} (W 1 ). 
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Let 



Jr"\b„ \x - y\ 



'R"\B H I a; - y| 

9R{x)= KJ ^dy, 

Jr"\b r \x - y\ 



then we know 



«(*) = / , ^La rfy + /«(*)■ 

JBr F - j/ 1 

«(*)= / , "^La ^ + gaW. 
■/Bh F - y| 

and /„ G L" +1 fl L,~ (B R ), g R G n ££ c 

To continue, we observe that by symmetry, we may assume p > q, then p > 
and p — ^ (p + 1) > 1. On the other hand, it follows from + ^j-j- + ^ = 1 that 
M-l = f(p+l)(g + l). Hence 

p--(p+ 1)1 5-l = -(p+l)>0, 
L n J n 

and this implies g^ 1 < p — ^ (p + 1). Choose r such that 

1 < r < p (p + 1) and < r, 

n 

for example, we may take r=p — ^ (p + 1), then 



v{*) 1/r <([ r^^dy) +9R{X ) 
v^-bk f y\ * 



1/r 



Ibr \x - y\ 

We have 



(*)=/ p(y }~ 

JBs F - 2/1 



/ f v(y) q ~ r / /" u(z) p - r u(z)\ \ ' u . . 

< / | ,n-a / ,n-V ^ + ftfl (*) . 

Js R kc - y V^Bn \y-z\ 



Here 



h R (z) = / : ^1=3 rfy + / fl (a:) . 

jb r F - y\ 

It follows from the fact that g R G (.Br) n (.Br) that h R G n 
i^ oc (Sij) for all q < oo. Let 

p+ 1 , q+ 1 



a = 



p — r q — r L 



then calculation shows — + £ = — + — , moreover we have 

ra b rn n ' 

rip, 
p + 1 a p + 1 

and 

1 _ a = P~g(p+l)-r ^ Q 
an p + 1 
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Hence for any p + 1 < q < oo, when R is small enough, it follows from Proposition 
12.11 fbv choosing a, j3, p, q, r, a, b, u, U, V and / in Proposition ^. II as a, a, p+X, q, 
r i fip j q^--i i u > uP ~ r , v q ~ r and h R respectively) that u G Lfl {B R /^j. Since every 
point may be viewed as a center, we see u G L* oc (M n ). This implies v G Lj£c (R n ) 
and then it G Lf£ c (R n ). Now observe that fR,gn G C°° (-Br), it follows from the 
usual bootstrap method that u,v G C°° (R™). The fact u,d G L°° (R n ) under the 
assumption u G L p+1 (R n ) follows from carefully going through the above argument 
and applying Holder's inequality when needed. Note that 

Xb 1 (x) XK-VBi (x) 
u = * v H * v . 



By interpolation we know v G L s (R n ) for all q + 1 < s < oo, hence G L s for 
2±i < s < oo. Since 2±1 < ™ it follows form the fact G L — ~ £ (R n ) and 

x '^„ B _ 1 l a:) G L — +£ (R n ) for e > small that it (x) -> as |x| -> oo. The fact 
i> (a;) — > as |x| — > oo follows similarly. □ 



3. All solutions are radial 

In this section, we will use the integral form of the method of moving plane 
(developed in [CLQ2] ) to prove the radial symmetry of solutions to the integral 
system. Such radial property was derived in [CLOU ICLQ21 [L] under the further 
assumptions that both p and q are at least 1. Our approach works for both this 
case and the case when p or q is strictly less than 1. We will need the following 
basic inequality: assume < 6 < 1, a > b > 0, c > 0, then 

{a + c) 8 -{b + cf <a e - b e . 

Indeed, for x > 0, let fix) = [a + x) — (b + xf , then for x > 0, f'(x) = 
9(a + x) 6 ' 1 -9(b + xf^ 1 < 0. The inequality follows. 
For £ G K m and s > 0, we denote 

(m 
i=i 

Proof of Theorem \l.l\ By Proposition ^. 21 we know u,v G C°° (R n ), u (x) — > and 
u (a;) — > as |x| — ► oo. It follows from Hardy-Littlewood-Sobolcv inequality that 
v G (R n ). Without losing of generality, we may assume p > q, then we know 
V > ancl P > Q • Hence we may find a number r such that 1 < r < p and 
q < r. 

For A G R, we denote H x = {x G R™ : Xi < A}. For x = {x u x') G R n , let 
x\ = (2A — xi, x'). We also denote u\ (x) = u (x\), v\ (x) = v (x\), 

Bl = {xEH x -ux (x) >u(x)}, 
B v x = {xGH x :v x (x)>v(x)}. 
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Note that by a change of variable, we have 



■(*)= / 



v(y) 



\x-y\ 

v(y) q 



— d y 



= [ V{V) na dy + ( 
J H\ \x - y\ n a Je 



v{y\) q 



Hence 



>h x \x - y\ 
v{y\) q 



Jh x \x - y\ Jh 



h x \x\ - y\ 
v{y) q 



—dy- 



y\' 

This implies 

u(x\) -u(x) 

= I (v(yx) q -v( y y)(- 1 — 

Jh x \\x-y\ 

In particular, for i e BJ, we have 
< u (x\) — u (x) 
< [ {v(y x ) q -v(yY) 



<h x \x x - y\ 



:dy- 



y\ n - a \xx-y\ r 



dy. 



1 



1 



in — a i ir 

x - y\ \x\ - y\ 



dy 



< 



< 



-dy 



qr v(y x ) q r (v {y x ) 1/r - v (y) 1/r ) \^dy. 

JBi v ' \x — y\ 



It follows from Hardy-Littlewood-Sobolev inequality that 

\u\ - u\ LP+1 ( B u} 

<c(n,a,q,r) v q x ~ r (vl/ r — v 1 ^^ 
< c(n,a,q,r) 



q-r 



L * (Bl) 



= c{n,a,q,r) Kl^+ 1( ^ 
On the other hand, for x & B x , we have 



1 / r l/r 

v x ' —V' 



L(«+i)r(BS) 

L(*+l)r(Bl) 



, x f u(yx) p J , f 

v ( x *) = / i 7^dy+ / 

JBi \x - y\ Jb 



u(y) p 



bi \x\-y\ 



-dy 



u (yx) p 



h x \bi \x-y\ 



) —dy + j ^ dy 

n—a i> I i in— a f 

Jh x - 



< 



+ 



bi \x - y\ 



u (y\) p , 

l^dy 



>x\Bi \x\-y\ 
u(yf 



Bi \x\ -y\ 



;dy 



Jh x \bi \x-y\ n ~ a JH x \Bi\x x -y\ n ~ a 
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Since 



v (x) 



JB" \x - y\ Jb 



Ibi \x - y\ 



bi \x\ - y\ 



U{vx) —dy 

n—a f 



Jh x \b? \x - y Jt 



u(y\f 



\ x ~v\ 

it follows that 

0<v(x x ) 1/r -v(x) 1/r 
u{y\f 



h x \bi \x\ - y\ r 



;dy, 



< 



L 



bi \x-y\ 



T , f u(y) p 



hi \x x - y\ 



I dy+ [ dy 

Jbi \x - y\ n - a Jbi \x x - y\ n - a % 

u(y x ) p/r u(y) p/r \ 



l/r 



l/r 



\x-y\ {n - a)/r ' \x x -y\ {n - a)/ \ 
u(y) p/r u(y x ) p/r 



\x-y\ { - n - a),T, \x x -v\ { - n - a),T i 



dy 



dy 



l/r 



l/r 



< 



< 2 



u(y\) p/r -u{y) p,r u(y) p ' r -u{y x ) 



sp/r ' 



i I \i n ~ a)/r ' 

V \ x -y\ 

[u x (yf r -u(yf r ) 

I in-a 

\x-y\ 



\x\ -y\ 

r . l/r 



(n—a)/r 



dy 



-dy 



< 



2p 



u\(y) p r (u\(y) -u(y)) r 



l/r 



dy 



\x-y\ 

It follows from Hardy-Littlewood-Sobolev inequality that 



l/r l/r 
v x ' —V' 



L(«+Dr(B3:) 



< 



2p 



f 

JB' 



u\(yf r (u\(y) -u(y)) r 



y\ r 



dy 



l/r 



L«+i(BS) 



< c(n,a,p,r) \u P x r (u\ — u) 1 



I l/r 

£+1 , . 
>L p (B«) 



p—r 



l/r 



\ L ^(B^ iUX - U)ril ^(^) 



<c(n,a,p,r) 
= c{n,a,p,r) |uaI2^i(b») K ~ u \l^(b 
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Hence we have 

K -«lip+i(B«) 

p — r 

<c(n,a,p,q,r) \u\\ L r p+1 ^ |«a|^+i(b«) K ~ u Il»+i(b») 



= c(n,a,p, g ,r)|n| L ; +l(2Aei _ B „ ) |t;| 9 L9+l(2Aei _ Bs) |« A - u^i^) 

p— r _ _ i 

< c(n,a,p,q,r) \u\ L r p+1 ( 2Xei _ B ^ ML+i^n) |«a - «li,+i( B «) • 

Here d = (1,0, ••• ,0). 

After these preparations, we will use the method of moving planes to prove the 
radial symmetry of the solutions. 

First, we have to show it is possible to start. Indeed, for A large enough, we 
know \ u \ L p+i( 2 \ ei -B u ) can ^ e arbitrary small, this implies that 

\u\ - u\ LP+1 ( B uj < ^ \u x - u\ LP+1 ( B ^ , 

and hence \u\ — u\ LP+1 r BU \ = 0. It follows that B\ = when A is large enough. 

Next we let A = inf {A G R : B u w = for all A' > A}. It follows from the fact 
u (x) — ► as | a; | — > oo and u (x) > for all x that A must be a finite number. It 
follows from the definition of Ao that u\ Q (x) < u (x) for x G H\ Q . We claim that 
u\ = u. If this is not the case, then since 



v(x)= / (u Xo {y) p - u {yf) 



v Xo (x) -v(x)= I (u Xo {y) p - u (y f) [ l -— - - 1 ) dy 



y\ \x\ -y\ 



and 



ua„ (x) -u{x)= j {v Xo (y) q - v (y) q ) (- l — - - ) dy, 



we see u\ (x) < u(x) for x G H\ . This implies x^Xex-B^ — ► a.e. as A t Ao. It 
follows that M iP +i(2Aei-B") ^ as A t A . Hence 

1 

\u\ - u\ LP+1 ^ < - \u\ - u\ LP+1 ^ 

when A is close to Ao- This implies = for A close to Ao and it contradicts 
with the choice of Ao- Hence when the moving plane process stops, we must have 
symmetry. Moreover, u\ (x) < u(x) for x G H\ when A > A . Indeed, for any 
A > A , we can not have u\ = u because otherwise u is periodic in the first 
direction and can not lie in L p+1 . Hence u\ < u in H\. 

By translation, we may assume u (0) = max ie R« u (x), then it follows that the 
moving plane process from any direction must stop at the origin, hence u must be 
radial symmetric and strictly decreasing in the radial direction. It follows from the 
equation that v has the same properties. □ 
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